Shape from sound: toward new tools for quantum gravity 
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To unify general relativity and quantum theory is hard in part because they are formulated in 
two very different mathematical languages, differential geometry and functional analysis. A natural 
candidate for bridging this language gap, at least in the case of the euclidean signature, is the 
discipline of spectral geometry. It aims at describing curved manifolds in terms of the spectra of their 
canonical differential operators. As an immediate benefit, this would offer a clean gauge- independent 
identification of the metric's degrees of freedom in terms of invariants that should be ready to 
quantize. However, spectral geometry is itself hard and has been plagued by ambiguities. Here, we 
regularize and break up spectral geometry into small finite-dimensional and therefore manageable 
steps. We constructively demonstrate that this strategy works at least in two dimensions. We can 
now calculate the shapes of 2-dimensional objects from their vibrational spectra. 

PACS numbers: 04.60.-m, 02.30.Zz 
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The much sought-after unifying theory of quantum 
gravity, see e.g., [1|, should also unify or at least very 
closely link the mathematical languages of quantum the- 
ory and gravity, namely functional analysis and differen- 
tial geometry. A natural candidate for this mathematical 
link is the discipline of spectral geometry. Spectral ge- 
ometry asks, for example, if or to what extent the shape 
of a thin object such as a vase can be described through 
its vibrational spectrum. For quantum gravity, the more 
general spectral geometric question is if or to what ex- 
tent the curvature of compact n-dimensional Riemannian 
manifolds without boundaries (which henceforth we will 
simply call 'manifolds') can be described in terms of the 
spectra of canonical differential operators on the mani- 
fold, such as the Laplace or Dirac operators. (Compact- 
ness here merely serves as an infrared cutoff that ensures 
that the spectra are discrete). Note that it is a different 
branch of spectral geometry that was popularized by M. 
Kac's paper Can one hear the shape of a drum? It 
asked to what extent the boundary of a flat membrane is 
determined by its spectrum, see e.g., 0-0]. 

If indeed at least euclidean spacetime curvature can be 
described in terms of the spectra of canonical operators, 
such as the Laplace or Dirac operators, this could be very 
useful for quantum gravity. This is because these spectra 
are invariants, i.e., they depend only on the Riemannian 
structure itself but not on gauge choices such as charts 
or frames. Therefore, the degrees of freedom of gravity 
would be cleanly identified as the relevant set of eigen- 
values. Their dynamics and quantization would be free 
of difficulties stemming from diffeomorphism invariance, 
see e.g., 0]. Even if the spectra uniquely determine the 
curvature only in a finite neighborhood of each manifold 
this could be very useful for the perturbative quantiza- 



tion of gravity. 

Concretely, spectral geometry studies maps, 5*, which 
map a set of curved manifolds into one or several spectra 
of canonical differential operators on the manifolds. Since 
any such S is highly nonlinear, the study of the invertabil- 
ity of S is hard. In order to make the calculations feasible, 
previous work has therefore generally involved stringent 
assumptions, usually assumptions of high symmetry of 
the manifolds in question. There have been intriguing 
results of spectra determining the curvature, as well as 
examples of isospectral non-isometric manifolds. For a 
review, see, e.g., Little is known, however, about the 
case that would be of most interest in physics, namely 
the generic case of manifolds without symmetries. 

In [1,01, it was proposed to make this nonlinear prob- 
lem tractable by breaking it up into an iteration of lin- 
earized problems. To this end, the idea is to consider the 
derivative DS of S. It maps small changes of the metric 
to the corresponding small changes of the spectra, a prob- 
lem which is solvable by perturbation theory. If DS is 
indeed invertible, one can then iterate small steps of cur- 
vature reconstruction to obtain for any given spectra the 
corresponding curved manifold(s). The existence of non- 
isometric but isospectral manifolds would then manifest 
itself in a path dependence of the manifold so obtained. 

Regarding the question of the existence of isospectral 
non-isometric manifolds, we remark, that, crucially, DS 
acts on changes to the metric that have generally to be 
assumed to come in scalar, vector and tensor types (as 
do, e.g., cosmological perturbations). This means that 
the spectra of the Laplacians on scalars, vectors and ten- 
sors will need to be considered so that DS can possibly 
be invertible. This fact had been overlooked before the 
paper and it explains why the Laplacian on scalars 
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was often found to be insufficient to capture the metric. 

Here, we will now address the key question if DS is 
indeed invertible in the generic case and if, therefore, 
the spectra uniquely determine the curvature in a finite 
neighborhood of each manifold. We begin with the in- 
tuition that DS can only be invertible if the number 
of degrees of freedom in perturbations of the metric at 
most matches the number of degrees of freedom in the 
spectra, except that both numbers are infinite, of course. 
To make the problem tractable, we propose to regular- 
ize the manifolds in question by approximating them 
through suitable geometric graphs embedded in a higher- 
dimensional manifold. This yields an S that maps from 
a finite-dimensional space of geometric graphs into the 
finite-dimensional space of spectra of Laplacians defined 
over the graphs. The basic problem therefore reduces 
to determining when the finite-dimensional matrix DS is 
invertible, while ensuring that the spectrum of the con- 
tinuum Laplacian(s) is recovered in the continuum limit. 

While this program for constructing DS~^ is relatively 
straightforward in principle, it is technically quite in- 
volved. This is because there are scalar, vector and ten- 
sor degrees of freedom in the perturbations of the metric. 
This means that our program generally involves Lapla- 
cians that act on the analogs of scalar, vector and tensor 
fields over geometric graphs that approximate manifolds. 
The theory of such Laplacians has not yet been fully de- 
veloped, see e.g., [loj . 

However, we can immediately carry out our program 
on 2-dimensional manifolds. Their metric possesses only 
scalar degrees of freedom and therefore also the met- 
ric perturbations of 2-dimensional manifolds can be de- 
scribed by scalar functions. Also, all scalar functions can 
be expanded in the eigenbasis of the Laplacian on scalar 
functions. Therefore, crucially, there is exactly one shape 
degree of freedom for each eigenvector of the Laplacian. 
This is as many as there are eigenvalues, counting degen- 
eracies. Notice that, for higher-dimensional manifolds, 
by the same argument, the number of shape degrees of 
freedom again cannot exceed the number of eigenvalues. 
We will now show that DS is indeed generally invertible 
and also that the infinitesimal steps can be iterated. This 
then allows us to calculate shape from sound. 

Construction of the geometric graphs. Our 
aim is to explicitly implement the above strategy for 2- 
dimensional manifolds. The first step is to regularize the 
problem to finite dimensions by approximating the man- 
ifolds by surfaces of polyhedra. To this end, assume that 
from the origin of a large number, A'', of half lines 
emanate in uniformly randomly chosen directions. Now 
let us assume further that N positive numbers {ri}f^i 
have been chosen. For each i = l...A^ we mark on the 
I'th half line the point, pi which is at the radial distance 
Ti from the origin. In this way, we obtain N points pi in 
]R? . Here, pi = rjCi and arc normalized direction vec- 
tors. Further assume a triangulation of the set of points 



{pi\ so that we obtain the surface of a polyhedron. Any 
star-shaped Riemannian manifold can be approximated 
to arbitrary precision by choosing A^ sufficiently large. 
Generalizations of this construction for non star-shaped 
manifolds should be possible similarly. 

We will work with arbitrarily high but finite precision. 
To this end, let us fix a large A^, a set of random directions 
of the TV half lines and a triangulation. We then allow 
only the A^ values {r^} to vary. In this way, we obtain 
an A^-parameter family of polyhedra that approximate 
star-shaped manifolds. The key question then becomes 
whether from the vibrational spectrum of these polyhe- 
dral surfaces one can tell their shape. 

In order to answer this question, we first note that the 
number of real degrees of freedom that describe the shape 
of a polyhedron in the family is the number of radii, r^, 
of which there are A. We need to compare this num- 
ber with the number, A^', of eigenvalues {Aj}j^]^ of the 
Laplacian of the surface of the polyhedron. While there 
are several definitions for Laplace operators on graphs, 
they all possess the same number, A^', of eigenvalues, 
namely N' = N. 

To see this, note that functions on the manifold are re- 
placed on the graph by functions on the N nodes {pi}fLii 
see, e.g., 
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The function space is therefore an A^- 
dimensional Hilbert space. The Laplace operator, L, on 
functions over the graph acts on T which means that it 
is a self-adjoint N x N matrix . The spectral theorem 
then implies that L possesses N' ~ N real eigenvalues. 

The geometric graph Laplacian. In order to study 
the information contained in the spectrum in detail, we 
need to choose an explicit graph Laplacian. The choice 
is not unique, much like there is no unique discretization 
of a first derivative. General requirements are that it 
must be self-adjoint and semi-positive definite [l2|. Here, 
for the triangular surface meshes that we work with, we 
also need locality, linear precision and positivity of the 
weights [l3| . Choices of discrete geometric Laplacians are 
discussed in We will use the standard discrete ge- 
ometric Laplacian for triangulations introduced by Des- 
brun et al. It obeys these conditions and is also 

widely used in computer graphics. 

To define it, recall that a graph, G, of a polyhedron 
consists of a set V of nodes or vertices, Vi, given by the 
points Pi £ M^, and a set of weighted edges E. By dy 
we denote the degree of the vertex v, i.e., the number of 
edges emanating from it. The weight matrix W has as 
matrix elements the weights Wij that are associated with 
the edge that joins the vertices i and j: 



Wi 



cot(7y ) -t- cot(/3y) 



(1) 



Here, jij and (3ij denote the two angles opposite the edge 
(i, j). The Laplacian is then defined through L — V ~W, 
where V = diag{vi,V2, . . . , with Vi = Wij. 
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Crucially, this graph Laplacian is also known to be ge- 
ometric in the sense that it converges to the ordinary 
Laplacian in the continuum limit. Concretely, the re- 
sults of [3l show that if a set of so-called normal poly- 
hedra converges in distance and in area then the geo- 
metric graph Laplacian converges to the Laplacian on 
the manifold with respect to the operator norm topology. 
The theorem applies because, for sufficiently large N, our 
star-shaped polyhedra will be normal graphs in the sense 
of [l6l |. meaning that every point on the polyhedra has 
a unique closest point on M, whose distance is smaller 
than a fixed upper bound (the so-called reach of M). 
This means that given a star-shaped compact Ricman- 
nian manifold M , our polyhedral approximations possess 
spectra that converge to the spectrum of the Laplacian 
on the manifold for TV — >■ oo. 

Now as is easy to verify, zero is always an eigenvalue 
of this graph Laplacian, with the constant vector be- 
ing the eigenfunction (as on the manifolds). Thus, only 
the eigenvalues {Ai}^2 can carry shape information (al- 
though if we allowed nontrivial topologies of the poly- 
hedra, the degeneracy of the eigenvalue 0, being the 
O'th Bctti number, would carry connectivity informa- 
tion). Also, the graph Laplacian is conformally invariant 
since it depends only on the angles between edges but 
not on their lengths, like the conformal Laplacian on the 
manifold. Therefore, without restricting generality, only 
the radii {?'i}^2 carry shape information, while ri fixes 
the overall size. The numbers of degrees of freedom of 
shape and sound appear therefore matched at — 1. 

Infinitesimal Inverse Spectral Geometry of 
Graphs. In order to determine if, therefore, small 
changes of sound indeed determine the corresponding 
small changes of shape, let us consider the nonlinear 
map S : M^-^ namely S : (r2,...,rAr) -> 

(A2, Aat). Its derivative, DS, with entries dXi/dvj 
is a (TV — 1) X {N — 1) matrix. det(D5') is a generic 
function and should therefore be nonzero almost every- 
where (in the sense of on a set of measure zero) , so that 
(DS*)"^ should exist. Almost all small changes of a poly- 
hedron's spectrum should uniquely determine a corre- 
sponding change of shape. 

Concretely, since the graph Laplacian is conformally 
invariant, we fix one point on the graph, ri. As discussed, 
this fixes the overall size of the polyhedron and leaves 
us with only A^ — 1 radii {fi}fl2 determine from the 
spectrum, which has A^ — 1 variable eigenvalues {Ai}f^2- 
Now suppose the shape is changed by slightly changing 
the radii {?'i}^2- 

n^r, + ea,, i^2,...,N (2) 
where e ^ 1 and all a,; G JR. Then, 

L — !> L + eA(a), with a = (02, ajv) (3) 
We find that A = 6V - 6W, where Svj = J^j Swij, with: 
5wij = 77^,j cot + ri^^^ cot (4) 



The indices 7 and /3 denote the number of that vertex 
opposite the edge («, j) which has the angle or Pij re- 
spectively. Here, T]ki;j = [kk ■ ftk + Ijk ■ fik)/{hk ■ Ijk) with 
hj = fiCi ~ rjCj and = UiCi — ajCj. The change in the 
spectrum of L can now be calculated by first order per- 
turbation theory. With fc^*) denoting the i'th eigenvector 
of L 

N 

A.->A, + e (5) 

i.e., the linear map DS reads 

{DS\j = d\/drj = (fc(')|A(a(^))|fc«) (6) 

where (a^-'-')fc = 6j,k is the perturbation vector in the 
j-direction. For any shape, i.e., for any point r = 
(r2, ...,rpf) in the (A^ — l)-dimensional space of radii, 3?, 
this yields DS{r). Its determinant, det(£' S'(r)) is a real- 
valued function over Jl. If there is no further loss of de- 
grees of freedom (such as that the zero eigenvalue never 
moves), the single equation det(DS'(r)) = should de- 
scribe a submanifold Z C 3? which is at most (A^ — 2)- 
dimcnsional. This then means that DS{r) is invertible 
almost everywhere, namely on'R — Z. Notice that spectra 
may uniquely determine shapes even on Z, similarly to 
how g : X ^ is invertible for all x, including x = 0, in 
spite of g'{0) ^ 0. 

To test that there is no further loss of degrees of free- 
dom (other than the conformal invariance and the con- 
stancy of the zero eigenvalue that we took care of) we 
should evaluate the determinant for all r G 3?. The ex- 
plicit form of det(I?S'(r)) is unwieldy due to its depen- 
dence on eigenvectors. It is straightforward, however, to 
evaluate det{DS{r)) numerically for any r G 3^. The 
results indicate that the determinant is indeed nonzero 
almost everywhere. This implies, by the inverse function 
theorem, that at least for almost all the polyhedra there 
is a finite neighborhood within which any change of shape 
is recoverable from its change of spectrum. 

Prom infinitesimal to finite inverse spectral ge- 
ometry. The fact that DS{r) is invertible almost every- 
where in the space of shapes, 3^, suggests to tackle finite 
inverse spectral geometry by iterating the infinitesimal 
steps. For example, a simple strategy for finding a tar- 
get polyhedron when given only its spectrum is to start 
with an arbitrary ansatz polyhedron. Then, one slightly 
perturbs its spectrum towards the target spectrum, while 
calculating the corresponding change of shape. The pro- 
cedure is iterated on a straight line path in the space of 
spectra, all the while calculating the change of shape un- 
til the target spectrum is reached. This works as long 
as det(DS'(r)) 7^ on this path, i.e., until Z needs to be 
crossed. In this case, one strategy is to temporarily re- 
duce the number of radii and eigenvalues that one varies. 
The resulting lower-dimensional submatrix of DS{r) can 
be chosen to have nonzero determinant and is therefore 
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invertible. Once Z is crossed, one can return to the 
straight path to the target in the space of spectra. 

Another strategy can be devised that is fundamentally 
more robust when Z needs to be crossed: We start again 
with an arbitrary ansatz graph, given by a set of radii 
{ri}^2- Fo"^ this algorithm, we use a distance function, 
(7, that measures the distance between this graph's spec- 
trum and the spectrum of the target graph. We perturb 
the graph, i.e., its radii, in the direction of the gradi- 
ent By iteration we obtain in the space of shapes, 
Ji, the steepest path to the target spectrum and shape, 
with respect to the chosen notion of distance. 

This algorithm does not automatically halt when cross- 
ing Z. To see this, for any graph, specified by radii 
r = {vi}, we define its spectral distance a to the tar- 
get graph through 
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be near isospectral, slowing down convergence. This is 
easily fixed by fixing several radii instead of just one. 

We have found that the algorithm allows one to recover 
graphs from their spectra all the way down to machine 
precision. In fact, in the generic case, i.e., when choos- 
ing random directions for the half lines, we have not 
found any iso-spectral non-identical graphs. Else, when 
the half-lines are chosen symmetrically, we do of course 
have that at least each graph and its mirror image possess 
the same spectrum. 




(7) 




where / is some arbitrary but fixed strictly monotonic 
function. More general a arc possible, of course, but this 
level of generality suffices. 

The only possibility for this algorithm to halt is when 
the gradient vanishes, i.e., when da/dri = for i = 
2, N. This happens when the global minimum of a{r) 



is reached at the target 



^target 



or at an isospectral 



graph if one exists. It can also happen at other critical 
points of a{r), where for all j = 2, N: 



da 



N 



- E2(/(A.M)-/(Ar'^*)) 



et,. Of (Kir)) d\{r) 



1=2 
N 



dXi drj 



drj 



N 



= E"'■(■D5')^J 



(8) 



Here, v^ = 2(/(A,(r)) - /(A,;*'"'^'^*)^^. Therefore, the 
algorithm can halt prematurely at a r e Z, but only 
if v{r) also happens to be in the kernel of [DS)^ . In 
this case, the algorithm can be unstuck by changing the 
definition of /, because this changes v{r). Intuitively, the 
algorithm flows down the landscape <j{r) over 31 along the 
gradient. It halts if it runs into a local minimum of (j{r). 
By changing / one can change the landscape so that the 
minimum and its catchment area are removed from r. 
Alternatively, one can keep / fixed and randomly move 
r sufficiently to escape the minimum's attractor. 

Numerical results. We implemented the gradient- 
based algorithm, choosing the half lines' directions ran- 
domly or regularly. We choose / : x — > l/x, so that 
the algorithm tends to first converge the smaller eigen- 
values, which determine the large scale shape, and then 
fixes the small scale structure. Technically, for large A, 
there exist polyhedra that are near translates and ro- 
tates of the target polyhedron and that therefore would 



FIG. 1: Snapshots of the algorithm starting with a sphere 
and finding the cube from its spectrum alone. 

Conclusions and Outlook. Methods that enable 
one to fully describe manifolds in terms of the spectra 
of their canonic operators should have great potential for 
quantum gravity. This is because these eigenvalues, being 
diffeomorphism-indcpcndent invariants of the manifold, 
can then be identified as the pure, gauge-independent 
degrees of freedom of gravity. If applied to spacelike 
hypcrsurfaccs, their dynamics may be quantized along 
the lines outlined for the Dirac operator in 0. Applied 
to 4-dimensional manifolds, euclidean quantum gravity 
can then be pursued by studying the thermodynamics of 
these spectra. In fact, the euclidean Einstein action in 
four dimensions is known to be expressible ver y s imply 
in terms of the eigenvalues of the Laplacian, @,ll3|- 

In 0, it was proposed, therefore, to break up the non- 
linear problem of spectral geometry into the iteration of 
infinitesimal linear steps. Here, we have made this idea 
workable, at least for two-dimensional star-shaped man- 
ifolds, by introducing a program of graph regularization. 
This showed analytically that small changes of spectrum 
gcnerically do determine the small changes of shape. The 
numerics even indicates that, in the absence of symmetry, 
the spectrum determines the shape globally. 

In order to extend our results to four (euclidean) di- 
mensions, it will be necessary to further develop the the- 
ory of Laplacians on tensors on geometric graph approxi- 
mations of higher-dimensional manifolds. But it may also 
be possible to pursue a similar program with ultraviolet 
regularizations other than graph regularizations, so as to 
preserve the smooth manifold structure throughout. 
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